function atanThetaFunc = atanThetaFunc(vec1,vec2)

    atanThetaFunc = atan2d(vec1(1)*vec2(2)-vec1(2)*vec2(1), ...
        vec1(1)*vec2(1)+vec1(2)*vec2(2));
    if atanThetaFunc < 0
        atanThetaFunc = atanThetaFunc + 360;
    end
end

